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ABSTRACT 

We demonstrate a new technique for detecting radio transients based on interferometric closure 
quantities. The technique uses the bispectrum, the product of visibihties around a closed-loop of 
baselines of an interferometer. The bispectrum is calibration independent, resistant to interference, 
and computationally efficient, so it can be built into correlators for real-time transient detection. Our 
technique could find celestial transients anywhere in the field of view and localize them to arcsecond 
precision. At the Karl G. Jansky Very Large Array (VLA), such a system would have a high survey 
speed and a bcr sensitivity of 38 mJy on 10 ms timescales with 1 GHz of bandwidth. The ability to 
localize dispersed millisecond pulses to arcsecond precision in large volumes of interferometer data has 
several unique science applications. Localizing individual pulses from Galactic pulsars will help find 
X-ray counterparts that define their physical properties, while finding host galaxies of extragalactic 
transients will measure the electron density of the intergalactic medium with a single dispersed pulse. 
Exoplanets and active stars have distinct millisecond variability that can be used to identify them and 
probe their magnetospheres. We use millisecond time scale visibilities from the Allen Telescope Array 
(ATA) and VLA to show that the bispectrum can detect dispersed pulses and reject local interference. 
The computational and data efficiency of the bispectrum will help find transients on a range of time 
scales with next-generation radio interferometers. 

Subject headings: Instrumentation: interferometers. Techniques: interferometric, pulsars: individual 
(Crab, B0329+54) 



1. INTRODUCTION 

Since the discovery of the first pulsar, radio obser- 
vations have steadily revealed a wide range of fast (< 
Is) transient phenomena. Pulsars (Hewish et al. 1968; 
Backer et al. 1982; Camilo et al. 2006), planets and ac- 
tive stars (Hallinan et al. 2007; GrieBmeier et al. 2007), 
and peculiar new transients (McLaughlin et al. 2006; 
Lorimer et al. 2007) all have distinct radio phenomenol- 
ogy related to the physics of the host and the medium 
through which the signal propagates. 

While large, single-dish telescopes have pioneered the 
study of fast radio transients, the nature of their design 
has several hmits. Single-dish telescopes can only localize 
an individual pulse to somewhere in its field of view (typ- 
ically arcminutes; Manchester et al. 2001). Timing of pe- 
riodic emission can localize very precisely, but individual 
pulses can't be localized well enough to easily identify 
optical or X-ray counterparts. Since interferometers si- 
multaneously have a large collecting area and a large 
field of view, they survey more efficiently than single- 
dish telescopes. Survey speed scales roughly the square 
of the product of the number of dishes and their diame- 
ter; this predicts that the VLA surveys almost 40 times 
faster than the Green Bank Telescope, assuming all other 
factors are equal. Despite this potential, next-generation 
radio interferometers (VLA, PAPER, LOFAR, ASKAP, 
and MeerKAT) have fast transients observing modes that 
are essentially identical to the single-dish concept: point 
one or more beams on the sky and search for pulses 
(Macquart 2011; Siemion et al. 2011). 

In principle, interferometers have much more informa- 
tion available to them because they detect the electro- 
magnetic wave over a distributed array of receivers. In- 



stead of forming a beam, interferometers can measure 
visibilities by correlating signals across the array. Visi- 
bilities contain spatial information that can reject inter- 
ference and localize a source to arcsecond precision over a 
large field of view. Traditionally, visibilities haven't been 
used in this way because it was difficult to run correlators 
this fast and the data volume was too large to handle or 
efficiently search for pulses (Law et al. 2011). As a result, 
a standard radio interferometer averages data to second 
time scales, losing nearly all information about millisec- 
ond variability in the sky. Even on these longer time 
scales, imaging algorithms are complex and subject to 
subtle statistical biases (Condon et al. 1998). The chal- 
lenge of data rates and volumes is expected to grow with 
next-generation interferometers, which will have more 
antennas, wider bandwidths, and larger fields of view 
(Per ley et al. 2011). 

Here, we present a technique that reduces the challenge 
of detecting transients in visibilities to the much simpler 
single-dish problem, while maintaining the utility of vis- 
ibilities. The technique is based on the "bispectrum", 
a quantity formed from visibilities that is sensitive to a 
point source anywhere in the field of view (Cornwell 1987; 
Kulkarni 1989; Rogers et al. 1995). Differencing visibil- 
ities makes it possible to interpret the bispectrum as a 
statistic for transients anywhere in the field of view on 
time scales from milliseconds to seconds. We describe an 
algorithm for real-time, arcsecond localization of tran- 
sients using raw visibilities from a radio interferometer. 

The theory of closure quantities and their application 
to transients is presented in §2. In §3 we test the theo- 
retical predictions with millisecond time scale visibilities 
from the VLA and the "Pocket Correlator" (PoCo) in- 
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strument at the ATA (Law et al. 2011). As shown in §4, 
this technique has a moderate computational demand 
and can reduce the flow of data to a manageable size. 
We conclude with thoughts about new science accessible 
with this technique and how it can address the growing 
challenge of big data in the study of transients with radio 
interferometers. 



2. THEORY OF TRANSIENT DETECTION WITH THE 
BISPECTRUM 

Closure quantities are the combination of visibilities 
from antennas that form a closed loop (Cornwell 1987). 
An example of a closure quantity is the closure phase 
or triple phase, which is the sum of the visibility phase 
from three baselines, 0ijk = (j^ij + (j^jk + 4^ki- This combi- 
nation is useful because it is independent of phase errors 
associated with any individual antenna. Closure quan- 
tities have largely been developed for very long baseline 
interferometry, because they make it possible to detect a 
source even if phase calibration is unstable (Rogers et al. 
1995). An interesting corollary is that the triple phase 
responds in the same way to a point source anywhere in 
the field of view. 

A related closure quantity is the product of visibilities 
on a closed loop, called the closure triple or bispectrum. 
The bispectrum can be written as bijk = aij ajk aki e^^^^j^, 
where aij are the amplitudes of visibilities for baseline ij. 
In the absence of signal, the bispectra have a mean at the 
origin of the real-imaginary plane with variance propor- 
tional to the cube of the thermal noise. A celestial point 
source is coherent across the array and moves all bispec- 
tra to positive real values proportional to the cube of 
its brightness (Kulkarni 1989; Rogers et al. 1995). This 
makes the bispectrum useful for the statistical detection 
of a source anywhere in the field of view. 

An interferometer can detect sources with a low sig- 
nal to noise ratio (SNR) per baseline, s, by coherently 
combining signals over many baselines. For a homoge- 
neous array of Ua elements, there are nu = na{na — l)/2 
independent baselines, so coherently combining all base- 
lines improves SNR to s y/nu- For combinations of three 
antennas, there are ntr = Tiaijia — l)(^a ~ 2)/6 closed 
possibilities. It can be shown that some triple phases are 
redundant (e.g., (/)i23 = ^124 + ^143 + (/>234), which leaves 
only (ria — 1) [jia ~ 2)/2 "unique" triples (Rogers et al. 
1995). However, the bispectrum is subject to different 
correlated noise properties, since it is formed by multi- 
plying complex visibilities rather than summing scalar 
phases. As a result, in the limit of small 5, all ntr triples 
are independent (Kulkarni 1989) and can be coherently 
combined to improve SNR. 

Interferometers can detect a source in many ways. 
Coherent beamforming sums complex visibilities with 
weights to direct the beam ^. Antenna-based incoher- 
ent beamforming sums the power detected at each an- 
tenna, while baseline-based beamforming sums power de- 
tected on each baseline (i.e., the amplitude of a visibil- 
ity; Thompson et al. 2001). In all of these cases, the 
SNR of a source detection can be written in terms of the 
SNR per baseline and size of the array (Kulkarni 1989; 
Rogers et al. 1995). The SNR for coherent beamforming, 

^ In this context, one can think of imaging as coherently beam- 
forming all possible beams. 



bispectrum, antenna-based incoherent beamforming, and 
baseline-based incoherent beamforming are: 
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Note that for constant sources, the relation for SNR5i5p 
is only valid when s is smaller than 1 (Kulkarni 1989). 
However, for transient sources the noise is measured off 
pulse. This removes the correlated noise effect and makes 
the scaling valid for all s. We refer to these relations as 
the "apparent SNR" , since they are valid when the statis- 
tics are measured relative to the standard deviation off 
source. For Gaussian statistics, the apparent SNR is re- 
lated to the false positive rate as ^ * ^1 — erf ^^^^^ 

(Rogers et al. 1995); a 5cr deviation corresponds to a 
false-positive rate of 2.87 x 10~^. 

An important complication to this SNR analysis is 
that an individual bispectrum measurement is not Gaus- 
sian distributed. Visibilities are Gaussian distributed, 
so multiplying three, noise-like complex visibilities to- 
gether produces a long-tailed distribution. The product 
of three independent rea/-valued random variables fol- 
lows the Meier G function Go'3 (Lomnicki 1967); we are 
not aware of an analytic expression for the real part of the 
product of three complex random variables. Regardless 
of the exact distribution, the Central Limit Theorem pre- 
dicts that the mean of many, independent bispectra will 
approach a Gaussian distribution (Rogers et al. 1995). 
In the presence of pure noise, each bispectrum in an ar- 
ray is independent (Kulkarni 1989), so the distribution 
of the mean bispectrum will approach a Gaussian as the 
number of elements in the array grows, as shown in Fig- 
ure 1. 

To include this effect in our sensitivity calculation, we 
simulated noise-like visibilities and directly measured the 
flux limit of the coherent beamforming and bispectrum 
statistics as a function of array size, n^. The flux limit is 
quantified by calculating s (from Equations 1) that would 
produce the simulated 5a false-positive rate under Gaus- 
sian statistics. This simulation shows that the flux limit 
of coherent beamforming matches the expectation for a 
Gaussian-distributed quantity. However, the mean bis- 
pectrum has a larger false-positive rate than expected for 
an apparent 5cr deviation, especially for small Ua where 
the Central Limit Theorem has less effect. We use these 
simulations to define an effective bispectrum 5a flux limit 
that produces the Gaussian false-positive rate. The ef- 
fective 5<j flux limit is higher than predicted by Gaus- 
sian statistics at a level of 53%, 26%, 13%, and 6% for 
Ua = 3,6, 12,24, respectively. For this work, we assume 
the mean bispectrum is Gaussian distributed, which is 
reasonable for arrays larger than a roughly a dozen el- 
ements for a 5<j flux limit. Future work will develop a 
more general expression for the false-positive rate and op- 
timizations for transient detection in the non-Gaussian 
domain. 

Figure 2 shows how the apparent SNR of four tech- 
niques scales with SNR per baseline for a 27-element ar- 
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Apparent SNR 

Fig. 1. — Simulated histograms of the noise distribution of the 
real part of the mean bispectrum. Each bispectrum is calculated by 
multiplying random complex visibilities for an array of na anten- 
nas producing n^^ independent visibilities. The real and imaginary 
parts of the visibilities are sampled from a Gaussian with standard 
deviation of 1 and mean of 0. The mean bispectrum for Ua = S 
(the smallest array with a bispectrum) is centrally peaked with long 
tails. For larger arrays, the distribution of the mean bispectrum 
narrows and approaches that of a Gaussian, as expected from the 
Central Limit Theorem. The distribution can also be used to esti- 
mate the false-positive rate for a signal of a given apparent SNR. 



ray like the VLA. The apparent bispectrum SNR scales 
as s^, giving it a strong response to bright sources and 
very little response to faint sources. However, the ef- 
fective bispectrum SNR, defined according to a Gaus- 
sian false-positive rate, would be less sensitive for large 
s. Our simulations suggest that the effective bispectrum 
SNR never exceeds that of coherent beamforming. 

The bottom panel of Figure 2 shows the effective 5<j 
flux limit for an array with Ua elements and sensitivity 
like the VLA. For the beamforming techniques, this plot 
essentially shows how the lines in the top panel intersect 
SNR = 5 as the array size changes. For the bispectrum, 
the effective flux limit is shown by correcting for sensi- 
tivity loss due to non-Gaussianity, as measured in our 
simulations. Generally, since large arrays are sensitive 
to pulses with SNR per baseline much less than one, co- 
herent beamforming and imaging become more sensitive 
than the bispectrum as increases. For large n^, the 
bispectrum and incoherent beamforming techniques scale 

similarly (s ex na ), so the bispectrum is always more 
sensitive than the incoherent techniques. 

To interpret the bispectrum as the brightness of a tran- 
sient requires subtracting the signature of constant emis- 
sion from the visibilities. Proper subtraction requires 
that the visibility fringe does not change over background 
subtraction time. The maximum fringe rate is equal to 
the product of the rate of change of the longest base- 
line in the {u,v) plane and the largest direction cosine. 
In units of radians per day, the maximum fringe rate is 

TT ^max ^f.o.v./A, whcrC 6rnax IS the loUgCSt bascliue, Of.o.v. 

is the field of view in radians, and A is the wavelength. 
More precisely, the error from subtracting visibilities 
leaves residual errors proportional to the brightness of 
sources in the field. Thompson et al. (2001) quantify the 
effect of time smearing as a reduction in the peak bright- 
ness of a source. For small changes, we can parametrize 
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Fig. 2. — Top: Apparent SNR for coherent beamforming, bis- 
pectrum, and incoherent beamforming techniques as a function of 
SNR per basehne for a 27-element array. For an effective SNR, 
defined according to the false-positive rate, the bispectrum SNR 
never exceeds that of coherent beamforming. Bottom: The effec- 
tive 5cr fiux hmits as a function of the number of antennas. The 
effective bispectrum fiux hmit includes a correction due to non- 
Gaussianity. The text marks the sizes of major new and planned 
observatories, including LOFAR (Stappers et al. 2011), MeerKAT 
(Booth et al. 2009), and ASKAP (Johnston et al. 2008). This as- 
sumes Stokes I detection with a SEFD of 400 Jy, 1 GHz bandwidth, 
and 10 ms integrations, which gives a la sensitivity of 63 mJy per 
baseline. For reference, the VLA 5a 10-ms sensitivity is 17 mJy 
for coherent beamforming, 38 mJy for the bispectrum, 50 mJy for 
baseline incoherent, and 85 mJy for antenna incoherent detection. 

this in terms of a 1.4 GHz VLA A configuration obser- 
vation: 



sp-sp 
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where r is the integration time scale. 

For visibihty (or image) subtraction, the peak flux lost 
from time smearing is similar to the flux residual in the 
differenced data given by Equation 2. This error will 
contribute to the mean bispectrum in proportion to the 
fraction of all triples that contain the long baseline. A 
single antenna is part of (n^ — 1) {ria — 2)/2 triples, so 
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errors for a long baseline are suppressed by the fraction 
of triples with that baseline, or a factor of 3/na. Table 
1 gives the longest time scale for which the bispectrum 
technique can be used on several radio interferometers. 
The time scale shows when the systematic errors intro- 
duced by background subtraction biases the mean bis- 
pectrum more than its thermal noise (Icr). This shows 
that most interferometers are suitable for using the bis- 
pectrum technique to search for transients as long as 1 
second. The VLBA is so extended that the fringe rate 
can only be subtracted on millisecond time scales^. 

The bispectrum also contains information about the 
spatial distribution of emission. For a point source with 
large s, the standard deviation of complex bispectra 
is as = VSs'^Q^^ where Q is the noise per baseline 
(Kulkarni 1989). Note that the standard deviation of 
bispectra is subject to correlated noise, as opposed to 
the standard deviation of the mean bispectrum in time. 
In contrast with a celestial source, interference is often in 
the near-field or subject to multipath propagation, which 
produces a random triple phase and a larger variance be- 
tween bispectra. Another way to think about this is that 
terrestrial interference can look like a spatially-extended 
transient. The variance between bispectra provides a 
simple way to quantify whether a transient is point-like, 
as expected. 

3. DEMONSTRATION 

To use the bispectrum for transient detection, we as- 
sume that the real part of the mean of all bispectra is 
related to the brightness of a single source in the field 
of view. To use the bispectrum to detect and localize 
pulses, we propose the following algorithm: 

• Subtract the mean visibility in time: Dif- 
ferencing visibilities removes all constant emission 
over some time scale, r. 

• Dedisperse visibilities: The frequency- 
dependent arrival time must be removed to 
maximize SNR per baseline before creating the 
bispectra. 

• Calculate mean and standard deviation of 
bispectra: The mean of bispectra gives a single 
value related to the significance of the transient, 
while the standard deviation is related to its spa- 
tial structure. Both the mean and standard devia- 
tion have dimensions of DM and time, so they can 
be used to set thresholds just as for a single-dish 
transients search. 

• Calibrate and image candidate pulses: Bis- 
pectra can be used to make an image, but cali- 
brated visibilities are needed for an accurate local- 
ization (Schwab & Cotton 1983). Visibilities need 
to be buffered long enough to capture the dispersive 
delay across the band. The pulse is then imaged in 
a separate process. 

Figure 3 shows how the bispectrum and coherent 
beamforming techniques detect pulses in PoCo data to- 
ward B0329+54. For both techniques, the plot shows the 

^ While this estimate uses the full field of view, the VLBA 
has a much smaller effective field of view due to the delay beam 
(Deller et al. 2011). 
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Fig. 3. — Top: Time series of the apparent SNR measured by 
the bispectrum (blue) and with coherent beamforming (red) in 5- 
antenna PoCo data of pulsar B0329+54. The PoCo data have a 
time resolution of 1.2 ms and have been dedispersed assuming the 
known DM of B0329+54. This panel shows calibrated data with 
B0329+54 at the phase center. Middle: The same plot but with 
the pulsar off the phase center. Bottom: The same plot but with 
no phase calibration. 

apparent pulse SNR when the pulsar is at the phase cen- 
ter, away from the phase center, and at the phase center 
but uncalibrated. The beamforming technique only de- 
tects the pulse when it is at the phase center and data are 
calibrated. The bispectrum detects the pulse regardless 
of its location or the phase calibration. 

The apparent SNR of the bispectrum and coherent 
beamforming techniques for data including several pulses 
are shown in Figure 4. The top panel in Figure 4 shows 
that the PoCo pulse SNR distribution measured with 
beamforming and the bispectrum follow the relation of 
Rogers et al. (1995) well. Since there are only nine base- 
lines and five triples in this analysis, the pulses shown 
here have relatively high s. The brightest pulse in this 
data has s ~ 3, corresponding to a brightness of about 
50 Jy. A similar search with a larger array would have a 
much slower rise in SNRbisp relative to SNRbfc- 

Figure 5 shows how the bispectrum can distinguish be- 
tween a celestial point transient and terrestrial interfer- 
ence. The top panel shows a typical dynamic spectrum 
for a VLA observation of the Crab pulsar with 12 mil- 
lisecond integrations. As is well known from single-dish 
observing, interference is bright enough to be detected 
even after dedispersion. The bottom panel shows the 
mean bispectrum and standard deviation over bispectra 
for the same pulse and interference. The interference has 
a large standard deviation of bispectra, while the true 
Crab pulse follows the expected theoretical relation for 
s > 1. 

4. DATA AND OPERATIONS RATES 

Another strength of the bispectrum technique is its 
computationally efficiency. Table 2 shows how the com- 
putational demand scales for three coherent fast tran- 
sient detection techniques at the VLA. These computa- 
tional demand values are scaled to operations per sec- 
ond (flops) by multiplying by the integration rate in Hz. 
For the bispectrum and imaging techniques, the visi- 
bilities are dedispersed for each baseline and DM trial. 
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Fig. 4. — Top: Apparent SNR of the bispectrum versus that 
of coherently beamforming of PoCo data toward B0329+54. The 
SNR is shown for each integration, regardless of significance. The 
dashed line shows the expected cubic dependence in the apparent 
SNR, as described in §2 and shown in Rogers et al. (1995). Bottom: 
A time series of the same two statistics for the same data. The 
rotation period of B0329+54 is about 600 integrations and four of 
the five pulses in this window are detected. For all detected pulses, 
the apparent bispectrum SNR is higher than that of beamforming, 
as expected for a small array as PoCo was used at the ATA. 

while coherent beamforming is traditionahy done be- 
fore dedispersion. The imaging FFT requires TVbi A^ch + 
5 TVpi^ log2Np-^ operations per image (Briggs et al. 1999). 
For the VLA D configuration, images require roughly 128 
pixels per side, so A^pix = 128. 

The bispectrum technique is the most computationally 
efficient technique for coherent pulse detection. Running 
this algorithm on 10 ms integrations at the VLA would 
require roughly 34 Gffops. Imaging is at least two times 
more demanding (this estimate excludes image cleaning) . 
The coherent beamforming is driven by the number of 
beams needed to fill the field of view. Even in the com- 
pact D configuration, roughly 7000 beams are needed. 
The VLA A configuration would require about 900 times 
more synthesized beams and 3 Pfiops to calculate. The 
large computational demand of coherent beamforming 
explains why it is typically only done with a few beams 
(Stappers et al. 2011). 

Our pulse detection algorithm consists of many inde- 
pendent processes (across frequency, DM trial, and time), 
so the problem is readily parallelized. Moreover, as- 




6 20 40 60 80" 



Fig. 5. — Top: VLA dynamic spectrum toward the Crab pulsar 
at 1.4 GHz with 12 ms integrations and 128 MHz bandwidth. The 
dynamic spectrum is formed from the sum of all cross correlations, 
effectively forming a synthesized beam at the location of the Crab 
pulsar with a brightness in units of Janskys. At the left is a pulse 
from the Crab pulsar, which appears as a red vertical line with 
a slight kink due to dispersion. Two channels show interference 
of similar strength toward the right. Bottom: The apparent bis- 
pectrum SNR versus standard deviation of all bispectra for each 
dedispersed integration. The blue dashed line shows the theoreti- 
cal expectation for a point source with s ^ 1. The one point near 
that line corresponds to the Crab pulse with s ^ 4. 

pects of the algorithm, such as the large number of mul- 
tiplications and dedispersion, are well suited to CPUs 
(Magro et al. 2011; Clark et al. 2011). Table 2 suggests 
that the bispectrum computational demand can be met 
with a single GPU-accelerated server. However, the data 
rate associated with millisecond integrations from a large 
array can limit its application for some correlator de- 
signs, so careful attention must be paid to data rate bot- 
tlenecks. For example, at the VLA, the visibility inte- 
gration time must be greater than roughly 10 ms, since 
nodes of the WIDAR compute cluster have 1 GbE input 
ports (Per ley et al. 2011). A correlator that can produce 
visibilities at high (> 10 Gb s~^) data rates is more likely 
to be hmited by computational resources, driven largely 
by the need to dedisperse visibilities from each baseline. 
In this case, it will be easier to run this algorithm in real 
time on arrays with a small number of elements. 

This description has focused on how the bispectrum 
can detect an astrophysical transient, but it can't be used 
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to localize it. While the bispectrum can be used to make 
images of the sky (Schwab & Cotton 1983), precise lo- 
calization requires external phase calibration. Thus, the 
final step in this algorithm is to save candidate transient 
data for calibration and imaging using traditional tech- 
niques (Briggs et al. 1999). This kind of post-processing 
is computationally more intensive than the bispectrum 
and is not feasible in real time on every integration. How- 
ever, since we only need to calibrate and image the occa- 
sional candidate event, the total computational demand 
will still be dominated by the steps before the event de- 
tection. 

5. CONCLUSION 

We demonstrate how to use the bispectrum to detect 
and localize transients with radio interferometers. The 
technique produces a single statistic proportional to the 
brightness of celestial transients, thus reducing the prob- 
lem of interferometric pulse detection to that of single- 
dish pulse detection. For a small cost in sensitivity, this 
technique makes an interferometer into a massive multi- 
beam receiver with a proportionately large survey speed. 
Using coherent beamforming to survey a similar area 
with next-generation radio interferometers would require 
10^ to 10^ times more computational power. 

The bispectrum technique is efficient enough to build 
into a correlator for real-time, millisecond transient de- 
tection. This system could save visibility data associ- 
ated with a candidate and apply standard calibration 
and imaging techniques to localize it to arcsecond preci- 
sion. The growing capacity of correlators for interferom- 
eters will make it possible to run such a system in real 
time to commensally search all data for transients. The 
technique can also be used to find transients offline and 
on longer time scales. In that case, the bispectrum is 
useful because it is computationally simple, resistant to 
interference, and requires no calibration. 

A real-time, commensal system would probe an un- 
precedented volume of data, giving it sensitivity to rare 
transients. Each pulse could be localized in real time 
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TABLE 1 
Time Scale for Using Bispectrum 
Technique on Radio Interferometers 



Telescope 


^f.o.v. 

(deg) 


^max 

(km) 


Error " 
(mJy) 


Tmax 

(s) 


VLA-A^ 


0.5 


36 


5 


0.2 


VLA-D^ 


0.5 


1 


5 


7.2 


VLBA^ 


0.5 


8611 


70 


0.003 


PAPER/ 


40 


0.3 


215 


1.6 


LOFAR^? 


2.4 


1 


40 


2 


MeerKAT^ 


1 


20 


7 


0.2 


ASKAP^ 


3 


2 


18 


1 



^Largest error expected in mean bispectrum, 
assuming Cassiopeia A located at the edge of 
the field of view (Baars et al. 1977). 

^Largest error expected in mean bispectrum, 
assuming Cassiopeia A located at the edge of 
the field of view (Baars et al. 1977). 

^Measured VLA sensitivity in A configuration 
at 1.4 GHz with 1 GHz of bandwidth. 

^Measured VLA sensitivity in D configuration 
at 1.4 GHz with 1 GHz of bandwidth. 

^Measured VLBA sensitivity using a 512 MHz 
bandwidth. 

■^Measured PAPER sensitivity at 150 MHz of 
100 mJy in 5 seconds with coherent beamform- 
ing of new 64-element array (Parsons et al. 
2010). 

^Predicted LOFAR sensitivity at 150 MHz of 
10 mJy in 1 second with 20-station coherent 
sum of core stations (Stappers et al. 2011). 

^Predicted MeerKAT "Phase 2" sensitivity 
with 80 dishes, 2 GHz bandwidth, and 30 K 
system temperature (Booth et al. 2009). 

^Predicted ASKAP "strawman" sensitivity 
with 30 elements, 300 MHz bandwidth, and 50 
K system temperature (Johnston et al. 2008). 



TABLE 2 

Computational Demand per Integration for VLA Fast Transient Search Algorithms^ 



Process 


Bispectrum 
Scaling Demand 


Coherent Beamforming 
Scaling Demand 


Imaging 
Scaling Demand 


Subtraction 
Beamform 
Dedisperse 

Image" 
Bispectrum 


lOA^biA^poiA^chA^t 14 Mop 
A^biA^poiA^chA^DMA^t 287 Mop 
16 A^trA^poiA^DMA^t 37 Mop 


lOA^biA^poiA^chA^t 14 Mop 
3ArbA^biA^poiA^chA^t 28 Gop 
A^bA^chA^DMA^t 2.8 Gop 


lOA^biA^poiA^chA^t 14 Mop 

A^biA^poiA^chA^DMA^t 287 Mop 
1.5e6A^DMA^t 602 Mop 


Total^ 


339 Mop 


33 Gop 


904 Mop 



"Excluding image cleaning, source identification, and application of calibration. 

^Assuming high-demand case: A^bi = 351, A^tr = 2925, Asb = 16, A^poi = 2, A^ch = 1024, A^b = 7000, and 
A^DM = 200. This produces spectra with 2 GHz of bandwidth and 2 MHz channel size. The number of dispersion 
trials, A^DM = 200, is appropriate for a high-demand case of DM = 1000 pc cm~^ for 1.2-2.0 GHz band and 10 
ms integration time. The number of time scales probed is effectively At = 2, assuming that the number of trials 
scales as the inverse of the time scale. Imaging and beamforming assume the D configuration, which has a longest 
baseline of 1 km. The A configuration requires 900 times more pixels and beams. 



